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Abstract 

Properties of a fundamental double-form of bi-degree {p,p) for p > 
are reviewed in order to establish a distributional framework for 
analysing equations of the form 

A^> + A^^- = 5 

(p) (p) (p) 

where A is the Hodge-de Rham operator on p— forms <1> on R"^. Par- 

(p) 

ticular attention is devoted to singular distributional solutions that 
arise when the source 5 is a singular p— form distribution. A con- 

structive approach to Dirac distributions on (moving) submanifolds 
embedded in is developed in terms of (Leray) forms generated by 
the geometry of the embedding. This framework offers a useful tool in 
electromagnetic modeling where the possibly time dependent sources 
of certain physical attributes, such as electric charge, electric current 
and polarization or magnetization, are concentrated on localized re- 
gions in spaced] 
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1 Introduction 



Mathematical distributions ^ play an all pervading role in the physical sci- 
ences. They underpin the theoretical formulation of quantum mechanics and 
many linear systems in classical field theory and are used to extend the so- 
lution space for linear partial differential equations to include non-smooth 
fields with singularities. They also offer a useful tool in electromagnetic 
modeling where the sources of certain physical attributes, such as electric 
charge, electric current, polarization or magnetization, are concentrated on 
localized regions in space [2], [3], [1], [5], [6] or spacetime. These aspects 
of distributions often occur together in problems in electromagnetic theory 
and this article explores this symbiosis in the unifying language of exterior 
double-forms. Such a formulation has direct relevance to many contemporary 
design problems in accelerator science where one is confronted with a number 
of difficult issues in electromagnetic computation. The effective resolution of 
such problems is often contingent on reliable computationally expensive nu- 
merical analysis which can greatly benefit from reliable analytic treatments. 
The formulation below will be illustrated by a number of applications in 
applied electro-magnetics. 

2 Notation 

The formulation exploits the geometric language of exterior differential forms 
[7] since this is ideally suited to accommodate local changes of coordinates 
that can be used to simplify the description of boundary and initial-value 
problems and naturally encapsulates intrinsic global properties of domains. 
In the following the distinction between smooth (C°°) forms on some regular 
domain and those with possible singularities or discontinuities will be impor- 
tant. Furthermore smooth forms with compact support will play a pivotal 
role. These will be referred to as test forms [8] , ^ , and distinguished 
below by a superposed hat. It will also be useful sometimes to notationally 

distinguish p— forms on different manifolds. Thus um will denote ap— formn 

(p) 

on some manifold M while 7 (X, Y) will denote a double-form of bi-degree 
{p, q) (see section [3] ) on the product manifold X x Y. Each n— dimensional 
^i.e. aM G SAp{M), the space of sections of the exterior bundle of p— forms over M 
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manifold M will be assumed orientable and endowed with a preferred n-form 
T*r 1 induced from a metric tensor qm- One then has [71 the linear Hodere 

(M) 

operator -k that maps p— forms to {n — p)-forms on M. If qm has signature 
tg one may write 

n 

9m = ^e'M0e^MVij (1) 
1=1 

where rjij = diag{±l, ±1, . . . ± 1) and 

^^^^ 1 = 4/ A A . . . A (2) 

with = detijjij) and {e\^} is a set of basis 1— forms in rT*(M). 
If a and /? are integrable p— forms on W C R" denote 



a^ -k (3 (3) 
by 

(a,/3)^^(Aa)^ (4) 

Furthermore if one of the forms has compact support on U (and belongs to 
the space of Schwartz test functions) an alternative notation is 

«^[0] = («,0)w (5) 
(p) (p) 

When U = R" the bracket subscript will be omitted. 

If (a, 0) is well defined (i.e. a is integrable with respect to the bracket (., .) 
) then one says that is a regular Schwartz p— form distribution associated 
with a. Not all distributions are regular and associated with smooth forms. 

Suppose Z is a smooth vector field on R" then Z = g{Z, —) is a smooth 
1— form. Using the property 

^^ = tgr-' (6) 

where the involution n on forms satisfies ria = (—lY a then it is straight- 

(p) (p) 

forward to verify the algebraic identities. 

^ ]=t,r/"-ia^[^-^ ] (7) 

W) {n-q) (.1) (n-q) 
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(ZAa^)[0] = a^M] (8) 

iza^[(3] = a^[Z Ap] (9) 

The operator S = d rj is a formal adjoint of the exterior derivative d 
with respect to the bracket (., .) in the sense that 



da[ ] = a[6(j)] (10) 
(p) (p+i) 



Similarly 



Sa[ ip ] = a[d^] (11) 
(p) (p-i) 

These relations are used to define da^ and 6a^. If / is a smooth 0— form 
then we denote by f the distribution defined by = <y^[f x]- 

If one defines the Lie derivative of with respect to a smooth vector 

(p) 

field Z by 

Cz OL^ = (iz d + diz) 
(p) (p) 

then 

£za''[x] = ayZAx + {Z'> + Z)6x-Vzx] (12) 
(p) (p) (p) (p) (p) (p) 

where V denotes covariant differentiation with respect to the Levi-Civita 
connection for which 6 = —g'^^ixa'^x^ 

In terms of d and 5 one defines □ the Hodge de Rham operator A = 
d5 + 5 d. A number of important identities follow by integrating the Leibnitz 
relation involving smooth forms 

d{a A P) = da A P + {r]a) A djS 

over regular domains U. These include 



(^0,7/-)^= (0,5V^)^+ / <P A^ij (13) 

JdU 



(50,^^)^ = (0,#)^- / ^ (14) 

JdU 



^The traditional Laplacian operator on forms is —A 
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Jau 

JdU 

(16) 

From these elementary identities further useful identities involving singular 
double-forms can be obtained. 

3 Double-Forms on 

The notion of a double-form was introduced in order to analyse problems in 
potential theory in pTj- Double- forms have since been used by Duff [12] 
and others as a powerful tool to study exterior equations on more general 
manifolds. A double-form of bi-degree {p, q) over the product manifold XxY 
may be regarded for each point in X as a p-form valued g— form on Y or for 
each point in F as a g-form valued p— form on X . To illustrate these notions 
suppose X = Y = R". In the natural coordinate system (a;*) for X and (?/*) 
for Y one hafl 

7(X,F)= E 7n...v,.i...ip(^'^)(^^^'A...Ac?x*^)0(rf|/^-^A...Af/y^'') 

(P''') il<...<ip jl<...<jq 

(17) 

The symmetric product is defined so that: 

{dx' dy^) A {dx^ dy') = {dx' dx^) A {dy^ dy') (18) 
and this implies 

7 A y = (-l)PP'+9'?' y A 7 (19) 
ip,g) {p',g') ip',q') ip,g) 

for all 7, 7'. The Euclidean metric tensors in natural coordinates and their 
associated Hodge maps will be designated: 

n 

gx = Ydx' (g)dx' (20) 

i=l 

Hi f : R" X R" R we denote its image by f(X,Y) with a similar notation for 
double-forms 
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gY = Y,dy'^dy' (21) 

i=l 

and 

* 1 = rfx^ A rfx^ A . . . A c/x" (22) 
(^) 

1 = c/yi A c/y^ A . . . A rfy" (23) 

The operators (ix, c^y, (^x, (^y, then act naturally on forms in the ap- 

propriate space. 

The fundamental double-form of bi-degree {p,p) is defined as 

7(X,F)= J2 f{X,Y){dx''A...Adx'^)xQ{dy''A...Ady'^)Y {24) 

^'P^ ii<i2...<ip 

where 

n>2 



{n-2)V„-i |x-y|"-2 

fiX,Y) = { (25) 
and@ 

K-i = K-2 / sin"-2^^rf^^ (26) 



^0 

with Vi = 271, Vo = 1 and |x — yp = Yl^=ii^^ ~ V^Y- This double-form is 
singular on X x F when x = y and satisfies 

7 (X,y)= 7 {y.X) (27) 
(p>p) (p.p) 

Ax7(^,>^)=0 X^r (28) 

dxl{X,Y)=5Y 7 (X,n (29) 
(p+i,p+i) 

5x 7 (X,y) = rfy 7 (X,F) (30) 
(p+i,p+i) 



is the volume of the unit rt-sphere in R" 



Ax 7 = Ay 7 (X,y) (31) 

(p.p) (p.p) 

The singularity at x = y in 7 (X, y) permeates many of the integrands 

(p,p) 

of integrals that feature below. Such integrals, when they exist, are then 
defined implicitly as limits. In general integrals a of n-forms a, over a 
domain U C R", that contain such a singularity are understood as 



Lim / a 

where is an n-ball of radius e centred on the location of the singularity. 
For any smooth p— form a on F the double-form 7 (X, Y) may be used 

(p.p) 

to define the type-preserving integral operator F: ax — ^ (r«)y mapping 
p— forms on X to p— forms on Y: where 

(F«)y= /«xA * 7 (^,>^)= / 7 (X,>^)A (32) 
7x (p,p) (^) 

From ([T3|) and f|T4l) it follows that for all test forms 0x, the differential 
operators and 6x commute with the integral operator F: 

dx{T^x) = Tidx^x) (33) 
Sxir^^x) = r(5x0x) (34) 



If one writes 



{T^)y=[ 7 {X,Y),<Px (35) 
(p) V(p,p) 



and notes that F is a symmetric operator: 

(Fa,/3) = {a,T(3) = (F/?, a) 

then it follows from Green's identity applied to forms with compact support 
that: 

= TAi! (36) 
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Thus 

(r0,Ax) = (0,r5x) = (0,x) (37) 

But(r(/),Ax) = (Ar0,x) hence 

(Ar0,x) = (0,x) (38) 

or 

Ar0[x] = 0[x] (39) 

and r is seen to be an inverse of A. 

For boundary value problems it is useful to introduce singular double- 
forms on lA that differ from the fundamental double-forms 7 (X, Y) by 

regular solutions to Laplace's equation A = 0. Such double-forms inherit 

ip) 

the singularity structure of 7 (X, Y) but can be tailored to satisfy certain 

(p.p) 

conditions on the boundary dU. 

As an example suppose one requires p— form solutions to the generalized 
Poisson equation 

A = 5 (40) 

that vanish on dU. Denote by 

g^^y^ 7 {X.Y)- ^^Y) (41) 
(p,p) {p,p) 

with 7 (X, Y) a symmetric non-singular solution of A 7 (X, F) = that 

(p,p) {p,p) 
coincides with the singular solution 7 {X, Y) on dU. Thus 

(p,p) 

Agp^y = X ^ r (42) 

^^y = ^^y (43) 

and 

g^xvlau = (44) 
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Suppose more generally now that Q^y is a double-form satisfying chosen 
conditions on dU. For any smooth p— form on R" define the integral operator 
Qhy 

{G<P)y = {G'xy.<Px) (45) 
If y G W the integral is singular at x = y and 

(GxY^'Px) = Lim {g^Y,<P)u-B, (46) 

One may use flTBl) together with the limits (derived by explicit computation 
in spherical polar coordinates with origin at the center of -B^): 



Lim / dxGxY ^ * 0x - 0x A -k dxQxY = '^^y 

Lim / I 6x(l)x A -k - SxY A ^ c^x^x I = 
where 5*^ is the n — 1 dimensional surface of the ball B^, to verify 

^4>Y = {Q ^4')y 

- j^^ {{H)x A ^^^G^xY - QIy a ^^{.dcfyx - SxG'xY A ^^^<Px + <Px A ^^G'xy^ 

(47) 

In this expression, = if y is outside W, = | if y G dU and u = 1 ii 
y G W If the first term on the right of this equation is written in terms of the 

source S in the equation (HOl) then ( 1471) offers a representation of the solution 

(p) 

of this equation in terms of this source and values of 50, d(j), -kcj) and cj) on 
dU. It should be stressed however that this representation is not an explicit 
solution in general since Scf), dcp, -kcj) and cannot in general be assigned 
values arbitrarily on dU. Nevertheless this representation is the cornerstone 
of many developments in potential theory [2]. The above identities and their 
subsequent uses are often attributed to Green. 

4 Electromagnetic Fields in Spacetime 

Maxwell's equations for an electromagnetic field in an arbitrary medium can 
be written 

rfF = and d ^ G = j, (48) 
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where F G TA^M is the Maxwell 2-form, G G Vh?M is the excitation 2-form 
and i E TA^M is the 3-form electric current sourc^. To close this system, 
"electromagnetic constitutive relations" relating G and j to F are necessary. 

The electric 4-current j describes both (mobile) electric charge and effec- 
tive (Ohmic) currents in a conducting medium. The electric field e G TA^M 
and magnetic induction field h G TA^M associated with F are defined with 
respect to an arbitrary unit future-pointing timelike 4-velocity vector field 
U G TTM by 

e = iuF and cq b = iu-kF. (49) 
Thus iuG = and iuh = and since g{U,U) = —1 

F = eAf7-*(cobAf7). (50) 

The field U may be used to describe an observer frame on spacetime and 
its integral curves model idealized observers. 

Likewise the displacement field d G TA^M and the magnetic field h G 
TA^M associated with G are defined with respect to U by 

d = iuG , and h/ cq = iu G . (51) 

Thus ^ ^ 

G = d Af/-^((h/co) At/), (52) 

and iud = and iuh = 0. 

5 Time dependent Maxwell Systems in Space 

In this article we restrict to fields on Minkowski spacetime which can be glob- 
ally foliated by 3-dimensional spacelike hyperplanes. The Minkowski metric 
induces a metric with Euclidean signature on each hyperplane. Furthermore 
each hyperplane contains events that are deemed simultaneous with respect 
to a clock attached to any integral curve of a future-pointing unit time-like 

^AU tensors in this article have dimensions constructed from the SI dimensions 
[M], [L], [T], [Q] where [Q] has the unit of the Coulomb in the MKS system. We adopt 
[g] ~ [i^], [G] = [j] = [Q], [F] = [Q]/[eo] where the permittivity of free space eo has the 
dimensions [Q^ T^M^^ L^^] and cq = ^ denotes the speed of light in vacuo. Note 

that, with \g\ = [L^], for forms a in n dimensions one has [ka] = [a][L"-^'^P] 
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vector field U = on spacetime and the Hodge map -k induces a Hodge 
map # on each hyperplane by the relation 

★1 = ccit A #1 

The spacetime Maxwell system can now be reduced to a family of exterior 
systems on R'^. Each member is an exterior system involving forms on 
depending parametrically on time t. Let the 3+1 split of the 4-current 
3-form with respect to the foliation be 

j = - J Ac/t + p#l, (53) 

(3) (2) (0) 



with ijiJ =0. Then, from (gHl) 

at (2) ^ 

dj = 0, (54) 

yields 

cZJ + p#l = 0. (55) 

(2) (0) 

Here and below an over-dot denotes (Lie) differentiation with respect to the 
parameter t and a = Ca_a. It is convenient to introduce on each hyperplane 
the (Hodge) dual forms: 

E := # e p := # d B := # b H := # h j := # J 

(2) (1) (2) (1) (2) (1) (2) (1) (1) (2) 

SO that the 3 + 1 split of the spacetime covariant Maxwell equations ( l48i) with 
respect to U = —c dt becomes 

de = -B, (56) 

(1) (2) 

dB = 0, (57) 

(2) 

dh = J + D, (58) 

(1) (2) (2) 

dp=p#l. (59) 

(2) (0) 
11 



All p-forms {p > 0) in these equations are independent of dt but have 

components that may depend parametrically on t. 

In the following it is assumed that b = u h and d = e e where e = e^en , 

(1) (1) (1) (1) 

it = UrIJ'Q ■ Thus in terms of e , h , E, H: 

(1) (1) (2) (2) 

d e = -ail, (60) 

(1) (2) 

dH = 0, (61) 

(2) 

<ih=eE+J, (62) 

(1) (2) (2) 

edE = p#l. (63) 

(2) (0) 

It is straightforward to show from this exterior system that the fields e and 

(1) 

h satisfy 
(1) 

A h + e/ih = 5J + # (^d e A #£-1 (^J - d hj^ (64) 

A e + 6/x e = -d (^Q - j + # (^d A # (^/."^ d ej ^ (65) 

Thus for homogeneous media with zero conductivity these equations each 
reduce to driven wave equations. Henceforth e, /i are assumed to be constant 

scalar s. 

Locally on spacetime one has an equivalence class of 1— forms whose ele- 
ments A differ by the addition of any exact 1— form. Then since F is closed 
one may write locally F — dA. Decomposing A (with respect to a particular 
U) into 0— form and 1-form A potentials yields 

A^-(l)U-cA (66) 

with iuA = 0. With U = —cdt and icl = #1 A U one has in terms of 
parameterised forms on 

e = -d(b - A 
(1) 
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b = #dA 

(1) 



The Maxwell system then reduces to 



5d(t) + 5k=^ (67) 
e 



5d A + e/iA = yU j — e/id (68) 



In a gauge with 

5 A - efKp = 

the Maxwell system above requires that the potentials must satisfy 

AA + e/iA = /ij (69) 



A(/. + e/i0=^ (70) 

for sources satisfying fl55|) . Physical solutions correspond to those satis- 
fying physically motivated boundary conditions (in both space and time). 
Alternatively in a gauge with 

5A = Q (71) 

the Maxwell system becomes 

= ^ (72) 



AA + e/iA = fi] — efi d 4> (73) 

In these equations on with the Euclidean metric, the Hodge map on 
all forms satisfies = 1 so 5 = ^d^r] and A = ^d^'qd + d ^d^r]. 
Henceforth, attention is devoted mainly to exterior systems on so it is 
unnecessary to maintain the notational distinction between d and d, S_ and 
6, A and A. 
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6 Electrostatics in with Sources in Do- 
mains with Boundaries 



Integral operators analogous to T may be used to solve electrostatic boundary 
value problems in involving the Hodge-de Rham operator A above and 
electric fields independent of t. As a basic example consider the Dirichlet 
problem of finding the electrostatic 0— form potential on C that 
must satisfy 

{o)_ (74) 
\du = <P 

for some data p = ^ EU,(j) G dU. 

The solution may be represented in terms of the integral operator Q de- 
fined by 



and is given by 



{Qp)y= / Px#QxY (75) 
Ju{x) (X) 

(t>Y = {Qp)y - I 0x # dxQxY (76) 

JdU (X) 



where 



gxY= 7 {X,Y)- 7 (X,r) (77) 

(0,0) (0,0) 

with 7 (X, Y) a symmetric non-singular solution of A 7 (X, Y) — Q that 

(0,0) (0,0) 

coincides with the singular solution 7 (X, Y) on dU. Thus 

(0,0) 

/^GxY = X^Y (78) 
QxY = Qyx (79) 

and 

QxY\du = (80) 
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If the data is smooth it defines re^uZar distributions p^, 0^ with supports 
on U and dU respectively by 

P^ix] = {P, X) ^""m = - f 0x # dx{gx)x (81) 

JdUiX) (X) 

Now 

{Gpm] = {Gp,x) = {p,Gx) = pm] (82) 

Thus for singular distributional sources one then has distributional 
Dirichlet solutions: 

<l>''[x]^p''m] + ^''[x] (83) 

If the eigenvalues Xm and Dirichlet eigenfunctions of A on can be 
found by solving 

(84) 

^M\dU — 

then a traditional way to satisfy the conditions required of Qxy is to express 
it as a Fourier expansion in the eigenmodes ^m{Y) for each point of F G W. 

If the real modes are labeled by a discrete index set M the expansion 
takes the form 

GxY = J2^^(^)^^(^) (85) 

M 

If such modes (with support on U) are ortho-normahsed so that ($m, ^n) — 
5mn then 

gN{X) = {gxY,^N{Y)) (86) 

Since the solutions $jv are zero on dU it follows that 
^n{Y) = Ajv / ^n{X) # GxY = AAr($iv(^), ^^xr) = XnQn{Y) (87) 

(X) 



or 



An 



Thus 



e?xy = J] -^^NiX) ^NiY) (88) 
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This series must be regarded as weakly convergent and the summations 
become integrations when the eigenvalues are continuous. 

For example consider U to be the space between two perfectly conducting 
parallel plates at z = and z = L, separated by a distance L in vacuo. Taking 
cartesian coordinates {x,y,z} for X and {x',y',z'} for Y the eigenvalues 
= {fc^, ky, n} with — oo < fc^ < oo, — oo < ky < oo, n = 1,2,... and (in 
complexified form): 

1/2 nnz , 

^N{x,y,z) = —\J — e-x.Y>[ik^ + iky) sm (89) 

This gives 
g{x,y,z] x\y\z') = 



sm sm 



47r2L^ L L J_^J^^ kl + kl + ^ 

(90) 

If the source p is not smoothly distributed in U one needs further technology 
to evaluate the solution and this will be developed below. 



7 Magnet ost at ics in with Smooth Sources 

Since magnetic charge is absent in Nature a typical problem in magnetostatics 
is the determination of a static magnetic field from a stationary electric 
current. The magnetostatic equations are a subset of Maxwell's equations 
and may be written: 

(5h=0 (91) 
(1) 

ci h = J (92) 

(1) (2) 

in terms of forms on R^. If the current is localized in space (vanishing at 
oo) the boundary conditions can be accommodated by using 7 (X, Y) for 

(1,1) 

the solution. If one sets |J h = ^dA./ jiQ then fPT]) is immediately satisfied. 



^Since rfB = one has B = c?A in reeular domains and in vacuo #B = unh where un 

(2) (2) (1) 

is the permeabihty of free space. 
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Writing J = — yUo#J one must have from (!92l) 
(1) 

(i#J = 

and in the gauge with 6 A = equation becomes 

AA = J 

Thus in free space with a smooth source jf 

A = rj 

since 

AA = AT J = J 



I.e. 



with 



Ay = /" 7 (x,r)A #:f 

Jx (1,1) (X) 

1 ^ 

(1,1) 47r|x-y|^ 



Note if i7 is bounded 



5yAy = [ 6y 1 {X,Y)Ai^Jx= [ dx 1 {X,Y) A#Jx 

Jx (1,1) (X) Jx (0,0) (X) 

= - / 7 (X,F)rfx#:^x = 

(0,0) (X) 



Since Ah = 63 one also has directly 

(1) 



h(F)= / 7 (X,F)A #(5J)x. 

.1) Jx (1,1) (X) 



(1) 

That this also furnishes a solution to (!9T|) and (p2|) will be verified in 

m. 



17 



Now for a smooth (regular) source J 

(1) Jy (y) 

= /"(/" 7 {X,Y)^ #jfx I A 
Jy \Jx (1,1) (x) J (y) 

= / / I 7 (X,F)A ] A 

JxJy (y) J (x) l-^uuj 

= / (rvi)xA #Jx 

7x (x) 

^ [ jxA# imx 

Jx (X) 

= j'^m 

If i7 is a smooth 1— form in the integral J^^^V"^] furnishes a solution 
for A. However if the current has support on a curve in it must be 
regarded as a distributional 1— form source not associated with a smooth 
1— form. Similarly the charge density in the electrostatic problem may be 
restricted to a curve or surface in R^ in which case a distributional source 

not associated with a smooth 0— form must be specified. 



8 Dirac Distributions on Submanifolds 

To define singular distributions with support on (possibly disjoint) subman- 
ifolds of a manifold M it is convenient to describe each submanifold para- 
metrically as an embedding in M. Suppose that a distribution has sup- 
port on a collection of submanifolds Sq \ i.e on the chain Ylir^o^- Recall 
that an n — A; dimensional submanifold S'g^'* (possibly with boundary) in an 
n— dimensional manifold M can also be prescribed in terms of k 0— forms 
/^i, /W2> • • • > f^"\ on M, such that rf/W^ A rf/^s • • • A df'^''\ ^ 0. Such 
forms generate a local foliation S^^'^ in the neighborhood of each submanifold 
Sq , where each leaf Sc of the fohation is an n — A; dimensional embedding 
given by f^'^\ — ci, f^'^\ — C2, . . . , f^'^\ — Ck for some constants ci, C2, . . . 
and c = {ci, C2, . . . Ck). We may choose the leaf Sq ^ with all these constants 
zero and denote it E'^^^^'' with = (/('■\, f^'^^, . . . , f^'^^). It is sufficient to 
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establish the distributional framework for a single component of a general 
chain so henceforth the label (r) will be omitted. 

If M is endowed with a metric g (and associated Hodge map the forms 
{dfj} give rise to a class of forms f2f on M defined with respect to -kl by 

(n-k) 

* 1 = d/i A d/2 A . . . A rf/fe A Vti (101) 

{n-k) 

This class restricts to a natural class of measures on Thus each 

representative induces the measure 



ujf = fif |^„-* (102) 

(n-k) (n-k) ^ 



on the submanifold Members of the class on M are equivalent if they 

differ by any combination of the dfj. This gauge freedom is of no significance 
on Ef "'^ since the dfj vanish there under pull-back. Thus it is sufficient to 
represent the class by an element satisfying 

tdf^ fif = j = l...k. (103) 

(n-k) 

In those situations where the leaves form orthogonal families, i.e. g~^{dfi, dfj) ■ 
for i j then 

-k {dfi A df2 A . . . A dfk) 

where = dfj (dfj) 7^ and the sign of \dfj \ is chosen so that the direction 
of the vector field dfj = g~^{dfj, — ) is in the direction of increasing values of 

■ df, 
-\dfj\- 



fj. Since dfj 7^ one has on E^ ^ two fields of "unit normals": Uj = ±t'^'^^ 



Then with n,- = jir- 



(n-k) \\dfj\ 



For non-orthogonal leaves an Qf can be chosen directly from fllOlj) . This 
condition is equivalent to the equation 

± 1 = i.{dfi A d/2 A . . . A d/fc A rif ) (106) 

{n—k) 
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which yields freedom to choose cuf with non-singular coordinate compo- 

(n-fe) 



nents when restricted to 



The singular Dirac 0— form with support on E^-'= is denoted dj.n-k and 

(0) 

defined by 



S^r.-k[(j)] = (j) UJf (107) 

(0) is;?-'= (n-k) 

Furthermore if F is a smooth vector field on M one may define the directional 
derivative of by: 

(0) (0) 

S^n-k = Y dj.n-k — iy dSj^n^k (108) 



(0) (0) (0) 



It follows that 



J) (0) ^sr' (n-k) 



5J„_40] = - / cuj (0) (109) 



where for orthogonal leaves: 

i^) = nfc "iwf|2 '^k ■ ■ ■ % d{i.$Y)\^n-k (110) 
(n-fe) ■'■■'■i=i 



Since 



6{f{Yij)) = f6{Yi;)-ij{Yf) 

(0) (0) 



(n~k) 



(111) 



Additional directional derivatives yield the distributions 

'(0) (0) 

For example 

Y^Y'^b^r.-k I [0] = (55.n-45(yi A 5(0^2))] (112) 

(0) ) (°) (0) 
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If one regards any point po G M as a 0— dimensional space then the 
original singular Dirac distribution 5^ with support on po is defined by: 

(0) 

S?M-^(Po) (113) 

(0) (0) 

It follows that 

<5E0=a;f<^S (114) 

(0) (0) (0) 

The above singular 0— form Dirac distributions can be used to construct 
singular p— form Dirac distributions with support on submanifolds in vari- 
ous ways. Thus in a spherical polar chart for a 1— form 0-layer Dirac 
distribution on the unit 2— sphere centred at the origin may be represented 

as 

{iie{eA)de + ii4>{e,4>)d4>) (115) 

(0) 

where iiq and /x^ are smooth functions on the sphere. A 1— form 1— layer 
Dirac distribution on the sphere (associated with some vector field Y on R^) 
may take the form 

0) de + p^{e, 0) (116) 

(o/ 

for smooth pg and p^. Similarly a 2— form Q— layer Dirac distribution on the 
sphere takes the form 

q{0,(f))d9 Ad(f) (117) 

m 

and a 2— form 1— layer Dirac distribution takes the form 

6^2s{e,(i))de Adcj) (118) 

(0)' 

for smooth q and s respectively. More generally in coordinates 

i^^, ■ ■ ■ ,^^,<^i, ■ ■ ■ ,<^p) 

n—p p 
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adapted to a p— dimensional submanifold in R", one has for j = 0, 1, . . .p, 
the J— form, "r-layer" singular Dirac distribution with support on Ef "'^ : 



(0) \ll<l2<-<l3 

(119) 

For a point electrostatic dipole with dipole moment Zx = P^-£j located 
at xo G X the electrostatic potential distribution 0y is given by the singular 
distributional source py where 

<Py[x] = Pm = i^z^ 7 (X, Y)] 5fjx] = P • V (t-^) x(xo) 

\ (0,0) J \|y-xo|/ 

(120) 



9 Line and Surface Currents in 

The source of a smooth magnetostatic field is the smooth current 2-form 
J in and the total electric current passing through a surface S C 'R^ 
is I[S] — Jg J, measured in amps in MKS units. In some circumstances 
this current density may be concentrated in the vicinity of 1-dimensional 
submanifolds Si C R'^ or 2-dimensional submanifolds S2 C R'^. Electric 
current filaments on material curves in space (wires) and current sheets in 
material surfaces in space are idealizations of such localized sources and may 
be modeled by singular distribution-valued 1-forms with support on curves 
and surfaces respectively. In these situations new physical current densities 
are introduced so that the total current in a segment of wire or a region of 
a surface is finite and measurable. These singular distributional sources can 
be mathematically modeled in terms of the distributions Ssr ior r — 1,2 and 
a vector field W in R^ with support on Ef . 

For a wire described by the curve with image Ej introduce the distribu- 
tional 1-form 

I^,^ = IoW (121) 
(1) ^ ' (0) 
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where Iq is a smooth function on the wire. Hence 

' (1) iR3(l) 



and with J" = # J: 

' ' Jsi (1) 7r3 Jr3 



El (1) Jsi (0) 

(122) 

since A cuf is zero on Sf . It does not make sense to ascribe a physical 
(1) 

dimension to since it is a functional not a value. However if {i7^['0]} 
denotes the physical dimension of J^[^] one has in MKS units 

U''[-]} = {^}{-} = amp{-} 

and {J} = {#J} =amp m-^ with {/o}{Vr}{c^f} = {J} = amp. The 

(1) 

(smooth) 0-form density 

i {s) = loiiwOJi) (123) 

(0) (1) ^ ' 

belonging to FA^Ej is the total current (amps) flowing along the wire at 
any point s G 

In a similar way, for the 1-form distribution on R^: 

^^^i^oW8^2 ^^24) 
(1) (0) 



one has 



(125) 



Ko{iwi')^f = / i> AiwiK'Qi^f) = / ip A k 
S2 (2) Js2 (2) Je2 (1) 



23 



since Ko^iw'ip) = '^oi' A (iw^f) on . Here 

(2) (2) 

(1) (2) 

is a smooth 1-form (measured in amps) on the surface Sj. If Ci G is 
any space-curve lying on this surface then j^^ k is the total current in amps 

crossing this curve in the direction W , on Sf. This may be compared with 
the definition above of the total current (in amps), J, crossing the surface 
S* G in the direction #J. If is a unit tangent vector |fl to Ci then 
ia_k is the surface current density in amp . 

It is worth noting that Jq G rA^R^, kq G ^A°R^ W G TTK? and Sj, E| 
are the primary objects used to define the distributions X^i and /C^2 ^ 

in terms of which are defined the smooth 0-forms / G FA^Sj and 1-forms 
k G TA^E^. 



10 Magnetostatic Fields from a Singular Dis- 
tributional Stationary Current Source 

Suppose a uniform electric current confined to a wire in R'^ is modeled in 
terms of a 1— form Q— layer Dirac distribution with support on the wire. 
A distributional source will generate a distributional 1— form A^. Thus 
define by 

A^[^]=jf^[r^] (126) 

Suppose J'^ has a support on the locus described by f = (/i, /2) with 
/i = r — a, and f2 = zm cylindrical coordinates (r, 0, z). This will be used 
to model a uniform stationary current source in a circular loop of radius a if 

= A(5si d(f) (127) 

(0)' 

^i.e. in terms of the Euclidean metric tensor g in R^, g{Ci-§^,Cl-^) = 1 
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with some constant A. Thus 

(0) 

^ ' (1) 



Y 



(128) 



Since 



and 



then 



(r^)y = / 7 (x,r)A (129) 

Jx (1,1) (^) 

= ^ (^.l(y) 7 ) A (130) 



A^[Vi]= / ^xA (131) 

(1) (X) 



where 

^x = A / Izs^(^) 7 ) (132) 

In cyhndrical polar coordinates, x = (r cos 0, r sin (f),z),y — (r' cos 0', r' sin 0', 2;') 
|x - y 1^ = i?2(r, 0, z] r', 0', ^') = + r'^ + z'^ + z'^ - 2zz' - 2rr' cos(0 - 0') 
and the Euchdean metric tensor is 

g ^ dr<^dr + r'^d(f)<^d(j) + dz<^dz (133) 

with 

_i d d _2 d d d d 
^ dr dr ^ d(f) d(f) dz dz 



Hence d(p' = ^-§^- With dx' = dr' cos 0' — r' sin 0' dy' = dr' sin 0' + 
r' cos (f)' d (f)' , id dx' = — r'sin^', is dy' = r'cos^', is d^' = 



Wy) 7 (^>^)lr'=a,^'=o = ^{c??^ sin(0-0') cos(0-0')} (135) 
-'(1,1) ait 
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Further with 



i^l = r' dr' A d z A d(j)' = d fi A d f2 A ujf 
(X) (1) 



cUf = adi 
(1) 



and one has |^ 



Ax{r,<P,z)=rd<p^ r \''''[^ ^] d<p' (136) 

Att Jq (r^ + + z'^ — 2ar cos(0 — (p'jj^'^ 

or with \1/ = 0' — and noting that the integrand is an even periodic function 
of ^: 

. / , N , , /"^^ Acos^d^ , 

Ax[r,(p,z) = -rdcf) — — ——j^ 137 

Jq (r^ + + 2^ — 2ar cos W)^'^ 

with magnitude independent of 0. For a constant current / in the circular 
loop the constant A = 47r /iq la. In MKS units the physical unit for J and 

hence h is the ampere and that for J is /io ampere/m. The methodology 

here illustrated for a circular planar coil is directly applicable to any open or 
closed current carrying conductor of arbitrary shape in space and possibly 
composed of piecewise smooth connected segments. 

To verify that when A-^ = J^A A jj^ip the distribution A-^ lies in the 
gauge satisfying 5A^ = one must compute 

SA^[i]=A''[di] 

(0) (0) 



(138) 



Y 



Aa 4^di = / di A4^A 

X (X) (0) Jx (0) (X) 



Jx (0) (X) Jx (0) (X) 

-A / e# / ^;^^^5x 7 {X,Y)u. 

Jxmx)JmY) "'^^^^ (1,1) (1) 



^The term in the integrand proportional to sin((/) — 0') integrates to zero. 
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But 

Sx 7 {X,Y) = dY 7 iX,Y) 

(1,1) (0,0) 

'^'^(^^ (0,0) d(f)' (0,0) 

so 

(5A^[e] = 

(0) 

since uj = adcp' and is a closed curve. This is physically equivalent 

(1) 

to the statement that the current in the coil is "conserved" . 

One can just as easily work in a gauge invariant manner by solving the 
equation 

Ah = 6J (139) 

All solutions to dH]) and ([22]) must satisfy ffT^ . As noted above for 
a smooth source 6 J, the solution in free space is 

h(Y)= [ 7 (X,F)A #(5J)x. (140) 

(1) Jx (1,1) (X) 

For a distributional source the solution 

fXoh''[i'] = J''ri'] (141) 
(1) (1) 

(cf f llOOp ) is modeled on the representation following directly from (11401) : 

f^oh^'li'] = [ / Jy A # 7' (X,F) A (142) 

(1) Jx Jy {Y){1,1) {X) 

In these equations the operator F' is defined as F but with 7 (A, Y) replaced 

(1,1) 

by 7' (A, Y) where 

(1,1) 

7'(A,F)=5y# 7 {X,Y). (143) 
(1,1) (^)(i,i) 

Thus for the distributional source 

(1) 
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for some vector field W on one has 

/io/^''[V^]= / Hx Ai^i'x (144) 

(1) (1) Jx (1) (X) 

where 

nx= [ XtWrl'{X,Y)ul (145) 
(1.1) (1) 

Following fll23p the choice of W and determines the physical dimensions 
of A. 

Although (11401) solves fll39p . in order to qualify as a magnetostatic field 
one must verify that it satisfies 6h = and dh = J. Now for some U G 
containing the smooth source J with 



hY={ 1 {X,Y),ix 

,(1,1) 



u 



one has 



{6h)y= {6y 7 {X,Y),ix =[dxl {X,Y) ,ix 



7 {X,Y),6xh + / 7 (X,r)A#xjx 

,(0,0) Jdu{o,o) 



iX,Y),ifxdJx] +[ ^iX,Y)AJx 
J ^ Jau (0,0) 

Hence with dJ = 0, and J =0 and (by explicit calculation) 



7 

.(0,0) 



Lim / 7 (X, r ) A Jx = 

(^^0) Jg^ (0,0) 

one has 

Sh = (146) 



-'^^For a loop with a non-constant but steady current, W ^ d4) and Yj\ ~ adip one could 
take A = 47r/io/oa/(^) for some physically dimensionless function / of to describe a coil 
with varying resistivity. 
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Writing = h and ■?/' = ■?/' in fll5l) 

But 5h = and ^5 J , -i/" j ~ ("^ ' '^^) ■ remains to calculate /^^ # rfh A 
with h = I 7 (X, r) , 5J I and so (rfh)y = ( 5x 7 (X, Y) , 5xJx ) • Thus 



L 



JidB.W \ (2,2) 



\ (2,2) 

where Ax = j,Q^ #y 7 (-^, ^) A -i/^y. In the limit as e — » one has by 

(2,2) 

explicit computation that this integral is zero, hence 

But since dip is an arbitrary test form it follows that 

dh = 3 on U (147) 
These arguments generalize immediately to the singular distributional case. 

11 Magnetostatic Fields from a Steady Heli- 
cal Line Current 

The derivation above of the potential for the magnetic field due to a circular 
line current was somewhat labored. However no further steps are required 
to find the potentials due to more complex geometries once one parametrises 
the source in terms of the geometry of its support in space. 

Suppose a helical wire with pitch p > and radius a is the space curve 

(r = a, = ^^^—^0", z = pa) in cylindrical polar coordinates. The parameter 
a is the arc-length parameter so < a < L for a helix of length L. It is 
convenient to introduce P > with a? P"^ + p^ = 1. Then the unit tangent 
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to the hehx is the direction W — + a P-§^- If one chooses to describe the 
hehx in terms of the fohating set f = {/i — r — a, f2 — z — ^ (j)} this yields a 
Leray form Q{ — —rd (f). Then, in this case, the source — XW dj^i with 

(0) 

support on the hehcal wire yields the magnetostatic potential 1— form 

A{r, 0, z) — Ar{r, (f),z)dr + A<p{r, (f),z)rd(j) + Az{r, (f),z)dz (148) 
where 

A = ^ ^ sin (Pa - 0) (149) 

= ^ cos (Pa - 0) (150) 



XapP f da 



Air 



with 



^ r"^ + + z'^ + a^ p'^ - 2 zpa - 2ar cos ( Per - 0) (152) 

For a constant current /q in the helix the constant X = Air Holoa. The po- 
tential A yields a global description of the field in terms of elliptic integrals 
that depend on the geometry of the helix specified by the parameters L, a 
and p or P. Such parameters offer natural scales that are useful for defin- 
ing dimensionless variables that in turn can be used to generate multipole 
or asymptotic expansions of the above integrals. The practical generation 
of high intensity uniform magnetic fields by carefully designing coils with 
complex helical windings is of paramount importance in the construction of 
undulators and free electron lasers. 



12 Magnetostatic Fields from a Steady Heli- 
cal Surface Current 

It is natural to model a solenoid composed of closely wound current- carrying 
turns by a surface current source. Such a current, regarded as a vector field 
on the solenoid surface Ef, can have an arbitrary direction H^|s2 and the 

distributions J = are well suited to model such a source. Suppose 

Ef is a right circular cylinder of radius a and length Lq and the surface 
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current is "painted" on it with helical strokes of pitch p = \J\ — P'^. Thus 
the integral curves oiW = p^+a P-^ are each similar to the helix above. To 
construct 6j^i one parametrises the solenoid surface as {r = a, z = p, (j) = cr} 
in cylindrical polars with < a < 27r and < p < Lq and takes f = r — a 
to ^nerate Qf = r d(f) A d z. It follows that, in this case, the source = 
Kq W with support on the cylindrical surface yields the magnetostatic 

(0)' 

potential 

A{r,(l),z) = Ar{r,(l), z) dr + Afj,{r,(l), z) r d(f) + Az{r,(l), z) d z (153) 
where 

.2n'^0 / / / "'P 



A^ = -a'P^ { ^cos{a-(f))] da (156) 
4vr ./n V./o R 



with 



p^ = jf- jf- + - 2zp - 2ar cos((t 



The double integrals above can be reduced to quadratures involving Elliptic 
integrals. 

13 Electrostatics with a Singular Distribu- 
tional Charge Source in a Domain with 
a Boundary 

Returning to the electrostatic example discussed in section [6] in Cartesian 
coordinates, suppose a uniform charged straight wire is inserted in the direc- 
tion of the y axis, between the grounded (0 = 0) planes at 2; = and z = L, 
at a position with z = zq {0 < zq < L) and a; = 0. Then with f = {x, z — zq} 
one has a 0— form 0-layer distributional source 

= A^si (157) 

(0) 
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and 



p^[gx]= [Vx#xx (158) 

Ju (X) 



with 



Vx = X QxY^f (159) 



El (1) 



Since ujj = dy' one finds from (M?]) 

(1) 



Vx {x, y,z) = — y sin — ^ sin ^ / ^e^^^^^ 160 

Tlie /cj: integration can be done by contour integration in tlie complex 
plane. One lias 




+62 6' 



xfe > 
xh <Q 



(161) 



yielding 



TT ^-^n n L L 

Vx{x,y,z)={ (162) 

Ay- isinnZLSOsinnp 2p q 
TT n L L 

This agrees with the computation in 



14 Time Dependent Electromagnetic fields and 
Smooth Sources 

When the smooth electromagnetic sources p and J depend on t they may 
generate time dependent electromagnetic fields that must satisfy (1^ . ([SI]), 
fl62|) . fl63|) . These equations offer a well-posed initial-boundary value problem 
that is traditionally approached by finding time-dependent potentials and 
A that satisfy the source-driven wave-equations fl69l) . fITOl) (subject to initial 
and boundary conditions) in the Lorentz gauge. When the sources are dis- 
tributional one seeks distributional potentials subject to similar conditions. 



32 



The distributional formulation given above for real static field configurations 
generalizes without difficulty to the time dependent situation. The essential 
modification is to define functionals on the space of t dependent complex 
test forms on and exploit the Fourier transform [13] of such forms. Thus 
if t/' is a t-dependent test p-form on X = its Fourier transform is the 
w-dependent test p-form J-'ip on X = R^ defined by 

{Ti)){X,u;) = j ij{X,t) exp{iujt)dt (163) 



^271 



R 



with inverse 



ip{X,t) = ^= / {J^4j){X,uj) exp{-iujt)duj (164) 



R 



The fundamental t-dependent {p, p) type double-form is the real part of the 
Fourier transform of Pi 

KUX,Y)=exp(-z-\^-y\) {X,Y). (165) 



It is useful, then, to define the singular complex-linear integral operator 
on all p— forms ip by: 

{T^ilj){Y,uj)= [ ir^(X,r) A #^(X,^) (166) 

Jx {p,p) (X) 

Suppose that jf is a smooth t— dependent p— form source (regular at spatial 
infinity) that enters into the equation 

AxCiX, t) + \ CiX, t) = JiX, t) (167) 

for the time-dependent p— form C on X = R^. Then 

(^C)(X,^) = (Ax - ^)-' {rj){X,uj) (168) 

"'^"'^Note that the phases in (|164p and (11651) are chosen so that for < r < oo, t > 
exp(±j^r). exp(=Fi(jji) — exp(±i^(r — ct)) describes a radially outgoing wave. 
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Since the singularity structure of K^^ (X, Y) is the same as that for 7 (X, Y) 

ip<p) (p,p) 
it is straightforward to show that a particular solution, regular at spatial 

infinity, is given by 

{TC){X,uj) = [ K^{X,Y) A # {TJ){Y,uj) (169) 

Jy {p,p) (Y) 

For each u the function JFC may be associated with a functional by applying 
the operation ^a{X,Lj) to define: 

(X) 

(J^C)[d]H= / (J^C)(X,cj) A #a(X,cu) (170) 
Jx (X) 

But from (I169p this may be expressed 

(j'CMiLu) = [ {rj){Y,uj) A # (r^d)(r,^) = (^j)^[r^d](^) (i7i) 

Jy (y) 
Thus from (11651) and Fourier inversion 

c{x, t)= [ 7 (X, F) A # :f f r, t - ^i^) (172) 

Jy (p,p) (Y) \ c y 



and 



c[a]{t) = ^ [ dt' [ :f(y,t')A# / 7 (x,r)A#(^d)fx,t-t'-^^') 

V27r Ji? (y) (p,p) (X) V c y 



(173) 



or 



c[a]{t) = ^ [ dt' [ :f('r,t-t'-^^V# / 7 (^,nA#(^«)(x,t') 

(174) 

Motivated by (11721) the distribution associated with the smooth p— form 
C (X, t) is defined by 

C^[P] = [ dt [ C(X,t) A # /3(X,t) (175) 
Jr Jx (x) 
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Hence 

C''0]=[dt[ [ 7 (X,F)A#jffr,t-^^)A#/3(X,y) (176) 
JR Jx Jy {p,p) (Y) \ c y 

or with t ^ t' — t — and interchange of X and Y 

= [ dt' [ {T^J){X,t') A # JiX,t') (177) 
Jij Jx (X) 

where 

(rK/3) (X, t')^ [ 7 (X, >^) A # /3 f F, t' + \^^^) 
Thus finally 

C^[/3] = / dt'J{X,t') A # (rR/3)(X,t') = (178) 

Ji?^ (X) 

where is the distribution associated with the smooth time-dependent 
p— form jf . If one now contemplates a distributional source J^^ that is not 
associated with a smooth time-dependent p— form on then (11781) offers a 
particular distributional solution to (11671) based on a retarded extension 
of the fundamental solution 7 (X, F). 

{p,p) 

The above solution is immediately applicable to the problem of finding the 
retarded electromagnetic potentials A and in free space, where e/z = 1, 
for any smooth time dependent source. Thus from section one sees that 
with a smooth source jf = — - in equation (I167P the smooth 0— form solution 
(j) associated with describes a free-space scalar potential solution to fITOl) 
and that with a smooth source S = A'-^J in equation (11671) the smooth 
1— form solution A associated with describes a free-space vector potential 
solution to fl^ . 

In the absence of losses due to_^ conduction the convective current space- 
time 3— form is given as j = po * ^ where the unit future-pointing time-like 
vector field V convects proper-charge density pQ. In a local spacetime co-basis 
{cdt,e^,e'^,e^} adapted to U the source velocity = 7 (^—cdt + ^) where 

7~^ = 1 — ^ and v = Ylk=i instantaneous time- dependent 

(Newtonian) 3— velocity 1— form in R^. Since 

J = -J Adt + p4^1 = Poi^V (179) 
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one finds immediately 

#J = -pv 

The analysis above also offers distributional solutions to fl69|) and flTOj) in 
terms of scalar and 1— form distributional sources. The electromagnetic field 
associated with an arbitrarily moving charge point source is space can now 
be modeled in terms of a moving Dirac distribution in space. 

15 Time Dependent Electromagnetic fields and 
Distributional Sources: Moving Dirac Dis- 
tributions 

The above section formulates a distributional solution to fll67p associated 
with a t— dependent p-form distributional source jf^ on R^. 0— layer sin- 
gular distributions associated with electromagnetic sources constrained to 
moving curves or surfaces in space can be constructed in terms of a family 
of distributions ^^^^^xj? defined for = 1, 2 by 

'(0) 



Then 




for some t— dependent 1— form W in with support on Ef 

For the case k = 3 the support of 6j.3-k is a moving point Po{t) in X = R^. 

(0) 

Then the worldline history of the moving point is encoded into the support 
of a distribution 6pgxR where 

(0) 

^oxR0] = [ Kpoit),t)dt (182) 

(0) JR 

If the history ol E X is parameterized by t it is the spacetime curve 

t'h^ (x,t) = (xo(t'),0 (183) 
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Then 



and 



Spoxr[^]= [ HMt'),t') dt' (184) 

(0) (0) JR (0) 

WSp,^R0]= [ {twit')P){Mt'),t') dt' (185) 
co^ (1) Jr roi 



(0) (1) Jr (0) 
One may now easily calculate 

Sp,M^rP] = [ {TRP){^o{t'),t')dt' 
(0) Jr 

^ I dt' I j{Mt'),y) #p(Y,t+ ^'^^^'^~y 

Jr Jy(o) (y) \ c 



(186) 



Changing variable t' t-^ t — t' + ^ at fixed Y with Jacobian Q defined 

by 

dt^ Q-\Y,t')dt' 

yields 

SpoM^rP]^ I dt f (Q(F,t07(xo(O,n ) 4hKY.t) (187) 

(0) JR JY \ (0) J {Y) 

where t' = f'(Y, t) solves the equation t' ^ t - ^^o^^^. Thus 

Z{Y,t) = Q{Y,t') 7(xo(t'),n (188) 

(0) 

is the 0— form associated with the distributional source Sp^^R [^rP] at the 

(0) 

field point Y at the instant t. It is straightforward to calculate the Jacobian 
in terms of xo(t) and its derivative xo(t): 

V c |xo(f)-y| J 

The instantaneous Newtonian 3— velocity of the point support is v{t) = Xo(i) 
and if one introduces the Euclidean unit vector 

n(F,t)- 



|y-xo(f) 
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connecting the field point y at time t to the source point :x.Q{t')a.t the earher 
time t' = t'{Y,t), then the Jacobian inverse takes the form 



Q-\Y,t)=(l-g(^,n{Y,t)] ] (190) 



where t' = i'{Y,t) solves the equation t' = t — 1^"*^*^)"^! and g denotes the 
Euclidean metric tensor on 

By contrast to the smooth electromagnetic sources discussed in section 
(|T4l) one models a moving point source (with electric charge q and 3— velocity 
v(t)) for the scalar potential (p^ by 

j'' = --Sp,y,R (191) 

(0) e (0) 
and for the potential 1— form by 

J^ = -g/iv(t)5p„xi?^ (192) 

(1) (0) 

It then follows immediately from (11781) that the distribution 0^ may be 
associated with 

V(Y,t) = -^-QiY,t) 7 {At'),Y) (193) 
(0) e (0,0) 

while the distribution A-^ may be associated with 

A(Y,t) = -g/iQ(y,t) M*') 7 W),Y) (194) 
(1) (1,1) 

where t' = i'(Y,t) as above. These are the classic Lienard-Weichert poten- 
tials for a moving point charge [2]. It is of interest to verify from i^J^ = 
—vp^ that the distributions A^ and 0^ satisfy the gauge condition S + 
efi(p^=0. It may also be noted that by using the appropriate fundamental 
double-forms in place of 7 and 7 these potentials maintain their struc- 

(0,0) (1,1) 

ture for solutions inVcY that satisfy appropriate boundary conditions on 
dD. 
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16 Conclusions 



With the aid of properties of the fundamental double- form of bi-degree (p, p) 
associated with the Hodge-de Rham operator A on differential forms, a dis- 
tributional framework for analysing equations of the form 

A$ + = S 

(p) (p) (p) 

on has been established. A set of r— layer Dirac singular distributions 
with supports on (moving) embeddings in R^ has been constructed and finds 
application in electromagnetic source modeling. The framework has been 
illustrated by explicitly calculating the fields associated with a current car- 
rying circular and helical coil, a finite length solenoid with a helical surface 
current density, a uniformly charged wire between two conducting plates and 
an arbitrarily moving point charge in free space. With the aid of r— layer 
Dirac singular distributions many fields with more complex distributional 
sources can be readily reduced to quadratures once one parametrises such 
sources in terms of the geometry of their support in space. 
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